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Abstract 


The scattering properties of a microstrip patch antenna with nonzero 
surface impedance are examined. The electric field integral equation 
for a current element on a grounded dielectric slab is developed for 
a rectangular geometry by using Galerkin’s technique with subdomain 
piecewise linear basis functions . The integral equation includes a resistive 
boundary condition on the surface of the patch. The incident field on 
the patch is expressed as a function of incidence angle. The resulting 
system of equations is then solved for the unknown current modes on 
the patch , and the radar cross section is calculated for a given scattering 
angle . Theoretical results in the form of radar cross section as a function 
of frequency are compared with results measured at the NASA Langley 
Research Center. 


Symbols 


L x 

dimension of the patch in 



x-direction 

d 

thickness of dielectric slab 

Ly 

dimension of the patch in 

pine 

^tan 

tangential components of incident 

y-direction 

electric field 

M 

number of subdivisions in 

pscat 

^tan 

tangential components of scattered 


x-direction 

electric field 

mn, pq 

indices specifying the subdomain 




basis functions on patch 


0 component of electric field 


number of subdivisions in 


N 

pjmn 

<P 

0 component of electric field 


y-direction 

j?inn 

Fourier transform of current 

R 

resistance matrix representing 


mode ran 


surface resistance on patch 

n 

dyadic Green’s function 

R s 

surface resistance on microstrip 
patch antenna 

jmn 

amplitude of mode mn 

T e 

characteristic equation for trans- 

J 

surface current on microstrip patch 


verse electric modes 


antenna 

T m 

characteristic equation for trans- 
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= v/=T 


verse magnetic modes 

K,a 

variables of integration in cylindri- 

yrrm 

component of excitation voltage 

cal coordinates 


vector 

K 0 

propagation constant for free space, 

{x,y, z) 

coordinates of field point 


2n j A o 

{x',y',z') 

coordinates of source point 

K x 

spectral domain transformation 

X 

unit vector in x-direction 


variable for x-direction 

x m j Vn 

coordinates of current mode mn 

Ky 

spectral domain transformation 
variable for y-direction 

Ax 

cell size in x-direction 
unit vector in y-direction 


y 

K x 

propagation constant for dielectric 
slab in 2 -direction 

Ay 

cell size in y-direction 


propagation constant for free space 

z 

impedance matrix to be solved 

K 2 



in ^-direction 

Zo 

impedance of free space, 377 £2 


2 

unit vector in ^-direction 

<r 

relative permittivity of dielectric 
slab 

O', <}>' 

incident angle of electromagnetic 
wave 

A 

piecewise linear function for current 
on patch 

Ao 

wavelength of electromagnetic field 
in free space 

O'O 

permeability of free space 

n 

pulse function for the current on 
patch 

°no 

#- polarized backscatter from 
^-polarized incident field 


0-polarized backscatter from 
^-polarized incident field 

a &0 

0- polarized backscatter from 
<p - polarized incident field 


</>-polarized backscatter from 
(^polarized incident field 

LJ 

radian frequency of electromagnetic 
field 


Introduction 

Spectral domain Green’s functions, which de- 
scribe the electric field radiated by a current source 
on a grounded dielectric slab, were introduced in 
the early 1980’s. This allowed the development of a 
moment method for analyzing perfectly conducting 
microstrip patch antennas. This technique accu- 
rately accounts for dielectric thickness, dielectric 
losses, and surface wave losses and can be extended 
to include the effects of a cover layer of a different 
dielectric constant on top of the antenna. Because 
of the spectral nature of the technique, it can easily 
be extended to model an infinite array of patches by 
examining only a single unit cell. Also of interest are 
the effects of lossy materials on the antenna. Lossy 
materials on the antenna will decrease the efficiency 
of the antenna, and hence the gain of the antenna will 
be lowered. This decrease in gain also means that 
the scattering from the antenna will be decreased. 
As losses are added to an antenna, other properties 
of the antenna, such as bandwidth, input impedance, 
and radiation patterns, will also be altered. 

The moment method technique incorporates ei- 
ther subdomain or entire domain expansion func- 
tions in order to model the current on the patch. 


Bailey and Deshpande (refs. 1 3) have used sub- 
domain expansion functions in order to model rec- 
tangular patches. Many other authors (refs. 4-13) 
have used entire domain expansion functions in or- 
der to model rectangular and circular patches. Bailey 
and Deshpande (ref. 14) have also used entire domain 
expansion functions to model an elliptical patch. 
The majority of this work has examined the input 
impedance and scattering properties of perfectly con- 
ducting patches both as single radiators and as in- 
finite arrays. Hansen and Janhsen (ref. 15) have 
included a space-varying surface impedance when 
modeling a microstrip feed network. 

A technique similar to the spectral domain 
method uses spatial domain Green’s functions with 
subdomain expansion functions in order to model 
inicrostrip structures. The disadvantage of this tech- 
nique is that it is not easily extended to examine in- 
finite arrays. A number of authors (refs. 16-19) have 
used this technique to model microstrip patch anten- 
nas as single radiators. Mosig (ref. 18) has mentioned 
that conductor losses can be included in this model, 
but no results have been presented. 

The boundary condition for the electric field on 
a thin resistive sheet has been examined by Senior 
(refs. 20-23) and is valid as long as the sheet is elec- 
trically thin. Using this type of boundary condition, 
several authors (refs. 24-27) have examined the scat- 
tering response of resistive strips and tapered resis- 
tive strips. This approach has also been used in order 
to study frequency selective surfaces (refs. 28 30). 
The same model for the surface resistance has been 
used in the study of superconducting materials and 
strip lines (refs. 31 and 32). 

This paper will describe spectral domain analysis 
of imperfectly conducting microstrip patch antenpas 
by using subdomain basis functions to model the 
patch current density. To simplify the analysis, the 
antenna feed will not be considered. The antenna 
is considered to be open circuited from the feed 
network, i.e., the feed impedance is infinite. Results 
are presented in the form of radar cross section as a 
function of frequency for a few representative cases 
and are compared with measured results. 


Theory 

The geometry of a rectangular microstrip patch 
antenna is shown in figure 1. The patch is assumed 
to be electrically thin and located on a grounded 
dielectric slab of infinite extent. The dielectric slab 
has relative permittivity e r , relative permeability /i r , 
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V 


where 




Figure 1. Geometry of a tnicrostrip antenna with arbitrary 
surface resistance. 


and thickness d. The standard time convention 
is assumed. The boundary condition on the patch is 
given by (ref. 23) 

EZi + - R*J (!) 


i r oc c to 

x e ji<s{x-x')ejKy{y-y') dKx dKy ( 4 ) 

and a and b can be x, y, or z. 

The components of the Green’s function are given 

by 


jZ„ A| A^K’jT’, + h^hyT^ 


Gxx (^fli Ky : d\(ij — K 


ifiTmT,, 


■ sin (K\ d) 


(5) 


~ , X )Z„ K * K V { K o T ”> - K \ K 2 T <) . IL . „ 

{K x .K v ,d\*) - V 7^. Sm(K|<i) 


(6) 


Gijx ( A".r, Ky, d\d) - G X y (A. r , d\d) (7) 


i/i/ (A^, A y . r/ </j 


j A’ 1 A 2 A y T,. -t A ,( A jT w 


sin (/v jrf) 


(8) 


G zx {Kt, Ky, d\d) = - { Z - A ..'. K 1 sin (A id) (9) 

rv J m 


The incident field is the field at the patch location 
due to an incident plane wave. The right side of 
equation (1) represents the field dissipated on the 
patch. The surface resistance R$ is, in general, 
a function of x and y and is equal to zero for a 
perfectly conducting patch. The scattered field is the 
field radiated by the electric current on the patch. 
Following a notation similar to that of Aberle and 
Pozar (ref. 9), the scattered field is found from the 
currents excited on the patch as 


G !y (K„K,Ad) = (ft, rf) (10) 

A q 1 m 

where 

T m = e T K 2 cos ( K\d ) + jK\ sin ( K\d ) (11) 

T v — K\ cos (K\d) + jK ‘2 sin (K\d) (12) 


gscat 




V ' 


z') • J ( c 7 , y, s') dx\ dy\ dz 

(2) 


where t? is the dyadic Green’s function for a current 
element on a grounded dielectric slab and J is the 
unknown vector electric current density on the patch. 
The dyadic Green’s function can be written in the 
form 


G — xGxxX + xGxyU + xG xz z 4- yG yx x 

+ yG vy y + yG yt 7 + zG zx x + zG zy y + zG 2Z z (3) 


K\ = \j (, Ko ~ J 2 (bn (A'! ) < 0) (13) 


K 2 = Jk* - P (In. (A- 2 ) < 0) (14) 


& = y/Kl + K* (15) 

The remaining terms of the Green s function are 
not needed in the present analysis. Details of the 
derivation of the Green’s function can be found in 
reference 3. Additional forms of the Green s function 
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are available in the literature and include such things 
as a dielectric covering above the antenna (ref. 10) 
and a uniaxial substrate (ref. 13). 

The current density J is modeled as a summa- 
tion of piecewise linear subdomain basis functions 
known as rooftop basis functions. This approach is in 
contrast to using entire domain basis functions that 
span the entire patch. Entire domain basis functions, 
such as sines and cosines, are useful for analyzing 
rectangular or circular patches but become cumber- 
some when used for other shapes. Mathematically, 
the subdomain basis functions for the components of 
the current density are described as 

M A r + 1 M N + 1 

■i. - E E J ""‘ = E E / ™ r ' A "* <*> ii " (v) ( is) 

1 71 = I m= l 7i= 1 


A/ + 1 JV A/+1 K 

J U = E E J y m EE ! »" A " <•"> <*> 07) 

If l — 1 77=1 77i=l 71=1 

where the functions A and n are “triangle” and 
'‘pulse” functions, respectively, and are expressed as 

{ 1 + (x - x m ) / Ax ( j'w - A x < x < x rn ) 

l (x ~ Tftt ) / Ax (x//i < x < Xjfi + Ax) 

M (Otherwise) 

(18) 


II" (V) 


1 

0 


(<)n - Aj/ < y < y n ) 
(Otherwise) 


(19) 


where Ax = 2 L X /(M -f 1) and Ay = 2 L y /(N + 1). 
Each single subdomain on the patch is specified by 
a pair of indices mn that also correspond to the 
coordinates x rn and y n . 

After equations (2) and (4) are combined, the 
order of integration may be changed and the basis 
functions that represent the patch current density 
may be transformed into the K x ,K y domain. These 
spectral domain current density functions are given 

by 


M ;V+1 

X (K x , Ky ) = E E W {Kxi Ky) (20) 

m= 1 71=1 


Af+1 N 

X (K x , K y ) j r F r (**. K y) ( 2 !) 


771=1 71 = 1 


where 


n"" (Ks. K y ) = Ax Ay 


sin (K y Ay/2) 

sin (A' x Ax/2) 

KyAy/2 

K x Ax/2 


x exp \-jK x x m - jK y y n + jK y ( Ay/2)] 


(22) 


( K , , K y ) = Ax Ay 


sin (K y Ay/2) 

2 

sin (K x Ax/2) 

KyAy/2 


K x Ax/2 


x exp \-jKjX m - jK y y„ + jK y (Ax/2)] 


(23) 


By using Galerkin’s method, the resulting equa- 
tions are then “tested” with the same set of basis 
functions, J p>l = Jx'x + J y q y, yielding a set of simul- 
taneous equations that can be solved by standard 
techniques. The boundary condition equation (1), 
becomes 


Jf JW • £tan dxdy=- II j pq . dx dy 

s s 

+ II J? q R s J mn dxdy (24) 

s 


and as p and q are varied over each subdomain, the 
resulting system of equations can be shown in matrix 
notation as 


i 1 

|?|? 

= 

~ j [gpqrnn-j ~ m 

r ypqmrn r ypqmn i 

[ *yx J L *yy J J 

1 1 


+ 

■ [ R p r n ] o • 
o [, Rj JT‘] . 

rmnl 

LL» J 


The impedance matrix terms are given by 
Z*T n = T1 [ [ Xx (k x , Ky,d\d) F” m (Kr.Ky) 

J-OG J OC 

X /f' ( h r. - Ky) <1K X dKy (26) 


z xy - . 2 f [ (k x , Ky,d\(i) F y in )K x ,K y ) 

J -OG J - DC 

X ( K x , ■ Ky) dK X dKy (27) 

Z yi = 712 [ j Gyx (k x , Ky,d\d) F" m [K x , Ky) 

J — OC J — OC 

X F” (-K x , -Ky) dK x dKy (28) 
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Z T" = T2 r r ( K *' K y'M) F ™ ( K ' K y) 

47r J — OC J - oo 

x F™ ( — K x , -ff y ) <*A'x ( 29 ) 

The integrations in equations (26)--(29) must be done 
numerically but can be simplified with the following 
change of variables 


K x = K cos q 
Ky = K sin a 


(30) 


With this change of variables, the integrals are 
changed to the form 

/ oc roc r2n roc 

/ [ ] dK x dK y = / / [ ] KdKda 

-oo J —oc Jo Jo 

(31) 

The integration from 0 to 2n may be further reduced 
to an integration from 0 to 7t/2 by using the even and 
odd properties of the integrand. It should be noted 
that each of the four submatrices of the impedance 
matrix are of a modified block Toeplitz form. This 
means that only the first row of each submatrix 
must be calculated by numerical integration. The 
remaining terms of each submatrix can be filled with 
the first row terms, thus reducing the time needed 
to calculate the impedance matrix. Furthermore, it 
can be shown that Zxy nn = Zy X ipq , which further 
reduces the computer time. 

The terms of the resistance matrix are given by 


ppqmn 

“xx ~ 


R» ( x , y) A p ( x ) 


rx p f Ax ry q 
J Xp — Ax Jy q — A y 

x II 9 (y) A m (x) n” (y) dy dx 


(32) 


jopqmn 

yy 


-r i 

Jx T) —Ax J ii 






R s ( x , y) II P (x) 


and equation (33) reduces to 




/? s Ax Ay 
— 

2Ks Ax Ay 



(q — n + 1 . q = n - 1 , p = rn) 

(q = u,p = rn) (35) 


(Otherwise) 


If the surface resistance R s varies across the patch as 
a function of x and y, the integrals may be evaluated 
by numerical integration. It is important to note that 
the terms of the resistance matrix do not depend on 
frequency. 

In order to examine the scattering from a micro- 
strip patch antenna it is necessary to evaluate the left 
side of equation (25). Each member of the excitation 
vector can be w'ritten as 


ypq = II '* Pq ' ^ int dx dy (36) 
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By reciprocity, this can be rewritten as (ref. 13) 


V1XI _ -4tt£W , Eg 

juHo 


(37) 


where E 0 is the vector amplitude of the incident plane 
wave, E pq is the far-field radiation from vector cur- 
rent mode pq on the patch, and —4n /ju?li 0 is the 
required strength of an infinitesimal dipole source to 
produce a unit amplitude plane wave. The factor 
e~jK° r J r has been suppressed in equation (37). The 
incident plane wave is assumed to have unit ampli- 
tude and is from the direction (0 l ,</> 1 ) in spherical 
coordinates with components Ep and E#. 

The fields radiated by a current mode on the 
patch can be found with the Green’s function defined 
above. The field at the point (x. y, z) from an impulse 
current source located at the point ( x f , y f , d) is given 
by 


x (x) A q (y) 11™ (x) A n ( y) dy dx (33) 


Note that the R p< x nU terms will equal zero if 
p > ( m + 1), p < (m — 1), or q ^ n. Likewise, the 
R pqmn terms will equal zero if q > (n + 1), 
q < (n — 1), or p ^ m. If R s is constant in the pq 
subdomain, equation (32) reduces to 


R™ mn : 


R , s Ax At; 

6 

2 R h Ax Ay 


{p = m + 1 , p = rn — 1 , q = n) 

(p — m,q — n) (34) 
(Otherwise) 


E a ( X , y, z) = ^ GatPK'W 

x e iKy{v-y') e -]K^-<0 dKx dKy (38) 

where a can be either x, y , or z and b is x or y . 
These integrals can be evaluated by the method of 
stationary phase and then integrated over the extent 
of each basis function to give the fields radiated by 
that basis function in the presence of the grounded 
dielectric slab (ref. 13). Once this has been done 
and the resulting equations converted to spherical 
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coordinates, the far-ficld components due to a single 
5-directed current mode, mode mn, are 


r -rmi 

b e 




v-jK 0 r 


(j K * d cos e 


x K\ A'„ cos <p p „ nl ^ ^ (3g) 


E"‘ n (r, 0. <j>) = 


Zo 

2-k 




°J K > d c.osO 


X A i S1I1 ^ S111 ( A l. (i l pnn {K ^ ( 4 0) 

T ( 


where K x and K y are evaluated at the stationary 
phase points: 


K x — — A" 0 sin# cos 0 
K y — —K 0 sin 0 sin <j> 


(41) 


Similarly, the fields radiated by a single ^-directed 
current mode, mode mn. arc given by 




~jh' 0 r ‘ 


e jKid c os (9 


K\ K 0 sin (f> sin (Kid) cmn t „ v , 

x (Ajc, Ky) 

■* 771 


(42) 


«*J' A * rf cos6> 

x A " ,os0sill,A|,/1 F” m (F x , A' v ) (43) 

where K x and are the same as in equation (39). 
By using equations (39)-(43) in equation (37), the 
left-hand side of equation (25) can be determined. 

Once the impedance matrix and the resistance 
matrix have been calculated, the results are added 
together to form a system of simultaneous equations. 
The excitation vector is then evaluated for the given 
incidence angle, and the system is solved for the 
unknown current coefficients. The scattered fields 
can then be calculated by summing the radiated 
fields from each mode on the patch at the given 
scattering angle. 

Results 

Computer programs have been written to evalu- 
ate the elements of the impedance and resistance ma- 
trices and then solve matrix equation (25). As men- 
tioned previously, only the first row of each Z pqnm 
matrix is calculated by numerical integration. The 


3T(r,0,*) = 


Zo 

2n 


-jh 0 r 


rest of each submatrix is then filled in by rearrang- 
ing the elements of the first row. Also, because the 
Z™T n and Zy X mn submatrices are related, only the 
submatrix is evaluated by numerical integra- 
tion. If results are needed over a band of frequencies, 
it is only necessary to compute the impedance matrix 
at a few widely spaced frequencies. The impedance 
matrix for other frequencies can be found by inter- 
polating each element of the impedance matri- 
ces. The number of frequency points at which the 
impedance matrix must be calculated by numerical 
integration depends on the span of frequency to be 
covered. This approach has been used with entire 
domain basis functions (ref. 12). As mentioned pre- 
viously, the terms of the resistance matrix do not 
depend on frequency and need to be calculated only 
once. The excitation vector, however, does depend 
on frequency and must be calculated for each fre- 
quency. After all the terms were evaluated for a 
given frequency, the system of equations was solved 
on the computer by using an IMSL library routine 
that solves a general complex system of equations 
and performs iterative improvement on the solution 
vector. 


When modeling the current distribution on the 
patch, it is necessary to chose M and N large enough 
to sufficiently approximate the true current distri- 
bution. Previous results (ref. 1) have shown that 
M = N — 7 is sufficiently large enough to model the 
current on the patch when the patch size is close 
to the size required for the first resonance. Exten- 
sive calculations have shown that for the swept fre- 
quency results presented, M = N = 12 is sufficient 
for modeling the patch current in the frequency band 
chosen. In figure 2 the calculated scattering from 
a rectangular microstrip patch is shown as a func- 
tion of frequency, where the scattering is measured 
in dBsm (dB/m 2 ). As indicated in the figure, there 
is no significant difference between M = N = 12 and 
M = N — 16. However, a slight difference is noted 
in the scattering response, depending on the num- 
ber of frequencies at which the impedance matrix 
was calculated by numerical integration. For the 
F = 3 responses shown, the impedance matrix was 
calculated at frequencies of 6.0 GHz, 10.0 GHz, and 
12.0 GHz. At frequencies between these three points, 
the impedance matrix was found by quadratic in- 
terpolation. The F = 5 responses were calculated 
essentially the same way with the addition of two 
more frequency points, 8.0 GHz and 12.0 GHz, where 
the impedance matrix was calculated by numerical 
integration. The addition of more points, i.e., F > 5, 
made no noticeable difference in the scattering 
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Figure 2. Calculated scattering from a perfectly conduc- 
ting rectangular microstrip patch antenna. L x - 0.75 cm; 
L y = 0.75 cm; d — 0.07874 cm; e r = 2.33; Loss tangent = 
0.001; (0\(^) - (60°. 180°). 


Figure 3. Calculated and measured scattering from a per- 
fectly conducting rectangular microstrip patch antenna. 
L s = 0.75 cm; L y = 0.75 cm; d = 0.07874 cm; e r = 2.33; 
Loss tangent = 0.001; (0\0*) = (60°, 180°). 


response. All the following calculated results have 
been performed with M — N = 12 and F = 5. 

Initially, the radar cross section of four microstrip 
patches, each with a constant resistance profile, was 
measured. The computed and measured responses 
for a perfectly conducting patch are shown in figure 3. 
The subdomain result was calculated as described 
above and agrees extremely well with the entire do- 
main result calculated by J. T. Aberle (Arizona State 
University, private communication). The result mea- 
sured in the Experimental Test Range (ETR) at the 
NASA Langley Research Center is slightly shifted in 
frequency and slightly lower than expected. This is 
not totally unexpected and can be attributed to the 
physical tolerances of the patch shape, dielectric con- 
stant, and dielectric thickness. The rapid fluctua- 
tions seen in the measured data are most likely due 
to imperfections in the background subtraction per- 
formed when processing the radar range data. These 
subtractions are necessary to approximate the re- 
sponse of the patch on an infinitely large grounded 
dielectric slab. Although the measured and calcu- 
lated data do not exhibit as close agreement as is 
evident in the subdomain and entire domain data, 


the relative position of the resonant peaks and the 
scattering levels are fairly close. 

The components of the scattered field for other 
polarizations arc shown in figure 4. The response 
shows a resonant peak at the same frequency as the 
ctqq response but also contains a peak in the center of 
the band where the oqq response does not. The upper 
resonant peak in the gqq response is not evident in the 
^00 response. The cross-polarized components, gq^ 
and <7^0, are both the same and show all the peaks of 
the previous two responses. The current density on 
the patch when illuminated with a ^-polarized plane 
wave as in figure 3, is shown in figure 5. The x 
and y components of the current are shown at the 
first and second resonances of the patch. At the 
first resonance the current resembles the expected 
sinusoidal distribution. At the second resonance this 
is also the case, although a whole period of the 
sinusoid is now evident. 

Calculations similar to those described above have 
been performed for a patch with a surface resistance 
of 5 over the entire patch. The calculated and 
measured results for <70# are shown in figure 6. Agree- 
ment in this case, is not as good as in the previous 
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Figure 

Lx 


4. Scattering from a perfectly conducting rectangular rnicrostrip patch as a function of frequency and polarization. 
= 0.75 cm; L y = 0.75 cm; d = 0.07874 cm; t T = 2.33; Loss tangent = 0.001; = (60°, 180°). 


6.25 GHz 


12.25 GHz 



case, but the general shape of the measured data is 
evident. Note that the peaks in the response have 
decreased and broadened compared with the per- 
fectly conducting patch. The complete set of scat- 
tering results are shown in figure 7. The decrease 
and broadening of the peaks is seen in each of the 
responses. The surface resistance was then increased 
to 11 n over the entire patch. The measured and 
calculated results for oqq are shown in figure 8. The 
agreement between the two in this case is quite good 
across most of the frequency range with some dis- 
agreement noted from 7.0 9.0 GHz. As expected, 
the resonant peaks in the response have decreased 


and broadened in shape. The current distribution on 
the patch is shown in figure 9. Although the gen- 
eral shape of the current distribution is the same as 
for the perfectly conducting patch, the amplitude has 
been considerably reduced. The a ee response for a 
patch with a constant surface resistivity of 20 ft is 
shown in figure 10. The calculated and measured re- 
sults for this case agree across most of the frequency 
band with only minor discrepancies at the lower fre- 
quencies. This is thought to have been caused by the 
measurement process, as evidenced by the rapid fluc- 
tuations in the measured data at the lower frequen- 
cies. With a 20 ft surface resistance on the patch, the 
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resonant peaks seen in the previous results are not 
evident. The radar cross section maintains a mono- 
tonic increase as frequency increases. The calculated 
results for gqq for these four cases are summarized 
in figure 11. As mentioned previously, as the sur- 
face resistance increases the sharp resonant peaks in 
the response gradually decrease and spread out. A 
patch with a surface resistance of 20 fl has no no- 
ticeable peaks in the scattering response. These re- 
sults suggest that the addition of surface resistance to 
the patch can be used to reduce the scattering from 
the patch and could possibly be used to increase the 
operating bandwidth of the antenna. However, this 
increase in bandwidth may be at the expense of low- 
ering the gain of the antenna. 


} ? 



Frequency, GHz 


Figure 6. Calculated and measured scattering from a rect- 
angular microstrip patch antenna with constant surface 
resistance of 5 ft. L x — 0.75 cm; L y = 0.75 cm; d = 
0.07874 cm; e r = 2.33; Loss tangent — 0.001; (6 l ,<p l ) — 
(60°, 180°). 

Additional calculations have been performed on 
patches with surface resistance that varies as a func- 
tion of position on the patch surface. A patch that is 
perfectly conducting but has a surface resistance of 
5 ft in each of the corners is shown in figure 12 along 
with the measured and calculated results for ctqq. 



Frequency, GHz 


Figure 7. Scattering from a rectangular microstrip patch 
with a constant surface resistance of 5 O as a func- 
tion of frequency and polarization. L x = 0.75 cm; Ly — 
0.75 cm; d — 0.07874 cm; e r = 2.33; Loss tangent = 0.001; 
(0\4> l ) = (60°, 180°). 



Frequency, GHz 


Figure 8. Calculated and measured scattering from a rect- 
angular microstrip patch antenna with constant surface 
resistance of 11 ft. L x = 0.75 cm; L y = 0.75 cm; d = 
0.07874 cm; e r = 2.33; Loss tangent = 0.001; (0\0*) = 
(60°, 180°). 


Close agreement between the two is seen across the 
entire frequency band, although the peaks in the 
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Frequency, GHz 
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Frequency, GHz 


Figure 10. Calculated and measured scattering from a rect- 
angular microstrip patch antenna with constant surface 
resistance of 20 H. L x = 0.75 cm; L y — 0.75 cm; d = 
0.07874 cm; e r = 2.33; Loss tangent = 0.001; {6 l , (p l ) = 
(60°, 180°). 


Figure 11. Scattering from a rectangular microstrip patch 
antenna as a function of frequency and surface resistance. 
L x — 0.75 cm; L y = 0.75 cm; d = 0.07874 cm; e r - 2.33; 
Loss tangent = 0.001; {0\4> l ) = (60°, 180°). 






6 7 8 9 10 11 12 13 * 


Figure 12. Calculated and measured scattering from a rect- 
angular microstrip patch antenna with surface resis- 
tance of R# = 5 12 for 0.375 cm < (|xj and |y|) < 0.75 cm. 
L x = 0.75 cm; L y = 0.75 cm; d = 0.07874 cm; = 2.33; 
Loss tangent = 0.001; (0* ,</>*) = (60°, 180°). 


measured response are not as high as was predicted. 
This discrepancy may have been caused by imper- 
fections in the shape of the patch. In order to en- 
sure that enough subdomains were used to model 
the current for this patch, additional calculations 
were performed with a higher number of subdomains, 
M = N — 16, and higher resistance on the patch cor- 
ners. These results are shown in figure 13 along with 
results for M — N — 12. As the resistance on the 
corners of the patch is increased, discontinuities in 
the current density may result. To accurately model 
this, a larger number of subdomains may be neces- 
sary. However, little difference in the results is seen 
in the figure with resistance in the corners as high 
as 100 12. The current distributions on the patch 
described in figure 12 at the two resonant peaks are 
shown in figure 14 and are similar to those shown pre- 
viously for patches with a constant surface resistance. 
A similar patch and the accompanying <tqq results are 
shown in figure 15. In this case the patch is perfectly 
conducting in the center and has a 5 f2 surface re- 
sistance around the perimeter. The resonant peaks 
in the response are lower than those of the previous 
response, as is expected with the addition of resis- 


Frequency, GHz 

Figure 13. Calculated results of scattering from a rect- 
angular microstrip patch antenna with surface re- 
sistance R s for 0.375 cm < (|x| and \y\) < 0.75 cm. 
L x = 0.75 cm; L y = 0.75 cm; d = 0.07874 cm; f r = 2.33; 
Loss tangent = 0.001: — (60°, 180°). 


tive material on a greater portion of the patch. A 
microstrip patch with a 5 f2 surface resistivity on the 
patch edges is shown in figure 16. The resistance in 
this case is on the patch edges that have the higher 
current density for the given excitation. A similar 
patch is shown in figure 17, but in this case the 5 Q 
surface resistance is on the patch edges that have the 
lower current density for the given excitation. It is in- 
teresting to note that the first resonance in figure 16 
is considerably lower than the first resonance in fig- 
ure 17, although the levels at the second resonance 
for each case are nearly the same. A similar patch is 
shown in figure 18, but the surface resistance on one 
of the patch edges has been increased to 20 fi. The 
predicted shape of the response can be seen in the 
measured data, although a frequency shift is clearly 
evident. As expected, with the increase in the sur- 
face resistance on the patch the level of the scattered 
field decreases. As a final example, a patch that has 
a 5 ft surface resistance above the diagonal and is 
perfectly conducting below the diagonal has been an- 
alyzed. This is shown in figure 19 along with the cal- 
culated and measured radar cross section a qq. The 
data show close agreement on the lower portion of 
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dBsn 



Figure 15. Calculated and measured scattering from a rect- 
angular microstrip patch antenna with surface resis- 
tance of R s = 5Q for 0.375 cm < (|x| or |y|) < 0.75 cm. 
L x = 0.75 cm; L y = 0.75 cm; d = 0.07874 cm; e r = 2.33; 
Loss tangent = 0.001; (0\ 4> l ) = (60°, 180°). 



Figure 16. Calculated and measured scattering from a rect- 
angular microstrip patch antenna with surface re- 
sistance of R s = 5 Q for 0.375 cm < \y\ < 0.75 cm. L x — 
0.75 cm; L y = 0.75 cm; d = 0.07874 cm; e r = 2.33; Loss 
tangent = 0.001; (0\4> l ) = (60°, 180°). 





Frequency, GHz Frequency, GHz 


Figure 17. Calculated and measured scattering from a rect- 
angular rnicrostrip patch antenna with surface resistance 
of R s = 5 Q for 0.375 cm < \x\ < 0.75 cm. L x = 0.75 cm; 
L y — 0.75 cm; d ~ 0.07874 cm; e r — 2.33; Loss tangent = 
0.001; {6\<P) = (60°, 180°). 

y 


Figure 18. Calculated and measured scattering from a rect- 
angular microstrip patch antenna with surface resis- 
tance of R s -50 for —0.75 cm < x < —0.375 cm and 
Rs = 20 H for 0.375 cm<r< 0.75 cm. L x = 0.75 cm; 
Ly = 0.75 cm; d = 0.07874 cm; f r = 2.33; Loss tangent = 
0.001; = (60°, 180°). 



Frequency, GHz 


Figure 19. Calculated and measured scattering from a rectangular rnicrostrip patch antenna with surface resistance of 
R S =5U above the patch diagonal. L* = 0.75 cm; L y = 0.75 cm; d = 0.07874 cm; c r = 2.33; Loss tangent = 0.001 ■ 
(0\ (p l ) = (60°, 180°). 
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the frequency band but have some differences on the 
upper portion of the frequency band. As before, this 
was thought to have been caused by imperfections of 
the patch shape. 

Conclusions 

This paper has described scattering from rectan- 
gular microstrip patches with a surface resistance Rs 
that is allowed to vary as a function of position on the 
patch surface. The boundary condition for the elec- 
tric field was used to derive an integral equation for 
the electric current that resides on the patch surface. 
Piecewise linear subdomain basis functions were used 
to model the current distribution on the patch sur- 
face. The necessary terms for representing the sur- 
face resistance on the patch were derived and were 
included in the equation in the form of a resistance 
matrix. The system of equations was then solved 
with standard solution techniques available on the 
computer. 

A variety of test cases were performed to ensure 
the validity of the theory and the accuracy of the 
computer codes. First, test cases were run for patches 
with a constant surface resistance. Additional cases 
were then run for patches with a surface resistance 
that varied as a function of position on the patch 
surface. Scattering results, ago, measured in the Ex- 
perimental Test Range at the Langley Research Cen- 
ter were compared with the predicted values. In 
all cases satisfactory agreement between the mea- 
sured data and the calculated data was noted. How- 
ever, slight differences were seen in some cases and 
have been attributed to physical imperfections of the 
patch shapes, the finite size of the ground plane used, 
and errors in the radar range background subtraction 
process. Scattering results for other polarizations 
have also been calculated and presented. 

The addition of resistance on the surface of a 
microstrip patch antenna has been shown to decrease 
the scattered energy from the antenna. A resis- 
tance of 20 ft on the entire surface of the patch to- 
tally removed the sharp resonant peaks evident in 
the frequency response of the perfectly conducting 
patch. In addition to reducing the scattering from a 
microstrip patch antenna, it may also be possible to 
increase the impedance bandwidth of the antenna 
with resistance on the antenna surface if the decrease 
in the antenna gain can be tolerated. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
May 24, 1993 
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